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Abstract 
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1 Introduction 



Charged black holes in four- dimensional N = 2 supergravity have attracted much attention in 
recent years. In particular, it is expected that the (indexed) entropy of = 2 supersymmetric 
black holes may be understood from a microscopic construction. From the macroscopic view- 
point, it is of much interest to understand the higher curvature corrections. Higher curvature 
corrections in regular Einstein gravity black hole solutions and their thermodynamic properties, 
have been studied as of the mid-1980's [T]. There, a solution in all space is obtained by treating 
the higher curvature term as a small perturbation. For higher curvature corrections of four- 
dimensional N = 2 supergravity, the horizon geometry and entropy of supersymmetric black 
holes have been studied a decade ago [21 [3] . The same have been studied recently, for extremal 
non-supersymmetric black holes [H [5] . Several works have considered the interpolation of super- 
symmetric horizon solutions in this theory to asymptotically flat space at infinity [HI [TJ El E] . In 
this work we consider N = 2 supergravity in four dimensions with small curvature correc- 
tions. We construct large charge extremal supersymmetric and non-supersymmetric black hole 
solutions in all space, and analyze their thermodynamic properties. 

Note that the i?^-terms considered in this paper are F-terms. One generally expects also 
D-teim corrections, which are not taken into account here. For supersymmetric black holes, it 
is conjectured that such terms do not contribute to the entropy [TO] . 

The paper is organized as follows: In section 2 we give a brief review of four- dimensional 
N = 2 supergravity with i?^-terms. In section 3 we discuss the framework and method for 
constructing the large charge black hole solutions. In section 3.1 we present the solutions and 
their thermodynamic properties in the case of supersymmetric black holes with i?^-terms. In 
section 3.2 we do the same for extremal non-supersymmetric black holes with i?^-terms. 

In the paper we will use a, 6, . . . = 0, 1, 2, 3 to denote the tangent space indices, corresponding 
to the indices /i, z/, . . . of the space-time coordinates (t, r, 0, 6). The exception are i,j = 1, 2 which 
are gauge SU{2) indices, and a = 1,2 which is a global SU{2) indices. The sign conventions for 
the curvature tensors follow [3]. 

2 i?^-Terms in = 2 Supergravity - A Brief Review 

We will consider N = 2 Poincare supergravity coupled to A'^v' Abelian N = 2 vector multiplets. 
N = 2 Poincare supergravity is a supersymmetric extension of Einstein-Maxwell gravity, adding 
two spin 3/2 gravitini to the graviton and (gravi-)photon. N = 2 Poincare supergravity can 
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be formulated as a gauge fixed version of = 2 conformal supergravity coupled to an = 2 
Abelian vector multiplet (see [3] for a comprehensive review). 

The on-shell field content of the vector multiplet is a complex scalar, a doublet of Weyl 
fermions, and a vector gauge field. We will consider A'^v' + 1 vector multiplets, and will denote 
by X^, 1 = 0.. .Ny, the scalars (moduli) in the vector multiplets. The couplings of the vec- 
tor multiplets are encoded in a prepotential F(X^), which is a homogenous of second degree 
holomorphic function. 

The N = 2 conformal supergravity multiplet (Weyl multiplet) is denoted by W"'^'^K It con- 
sists of gauge fields for the local symmetries: translations (P), Lorentz transformations (M), 
dilatations (D), special conformal transformations (i^), U{1) transformations (A), SU{2) trans- 
formations iy), and supertransformations {Q,S). In the theory without i?^-terms, the Weyl 
multiplet appears in the Lagrangian through the superconformal covariantizations. In order to 
get the i?^-terms, one adds explicit couplings to the Weyl multiplet. This appears in the form of a 
background chiral multiplet, which is equal to the square of the Weyl multiplet W"^. The lowest 
component of the chiral multiplet is a complex scalar denoted A. The prepotential F{X^ , A) 
describes the coupling of the vector multiplets and the chiral multiplet. One introduces the 
notation: 

^ ^F(XM), F^^ ^F(XM) , (1) 
and similarly for higher order and mixed derivatives. 

The bosonic part of the N = 2 conformal supergravity Lagrangian is 

Sne-^C = -i(X^Fi - X^Fi)(-R- D) + 

6 

+ {iV^X'V^Fj + '-FuiF-,^ - \xX-){F^'-^ - ^X-^T^'-) + 

- \FaA%b^' - '^K.F"'-) + -^f^^'-FaAK,' - \^'f:,) + 

- \%,F^jY^^' + '-F^C + h.c.) + L,^ . 

(2) 

e = A/|det(5'^,y)| where g^^v is the curved metric, K is the Ricci scalar, D is an auxiliary real scalar 
field of the Weyl multiplet, is the covariant derivative with respect to all superconformal 
transformations, V>a is the covariant derivative with respect to P, M, A, l^-transformations, 
F~^ is the anti-selfdual part of the vector field strength, is an anti-selfdual antisymmetric 
auxiliary field of the Weyl multiplet, and Yl, are real SU{2) triplets of auxiliary scalars of the 
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vector multiplet. The hatted fields are components of the chiral multiplet W"^, with their bosonic 
parts given by 



^0 A = T~^T^^- 



9"^ Bij = —lQR{V)(ij)abT"'^ 
F"^- = -lQn{M)J^T^- 

C = QAniM)j'-n{My''-, + 32R{V)jTR(VT'''; -IGT'^'-DaD'Tj, 



(3) 



= is the selfdual counterpart of the auxihary field, RiV)^fJ'j is the field strength of 
the SU{2) transformations, TZ{M)^i^^'^ is the modified Riemann curvature and TZ{M)^^'^~ is the 
anti-selfdual projection in both pairs of indices. The bosonic part of 7?,(M)^j'^^ is given bjll] 

= Ra,'' - 44^^^; + ^{T-.T'^'^ + TXT^'-) , (4) 

where -R^^^'^ is the Riemann tensor, and f^'^ is the connection of the special conformal transfor- 
mations, determined by the conformal supergravity conventional constraints, with the bosonic 
part^: 

fa' = - \{D + \RWa + \'R{A): + ^T;,T^'^ , (5) 

where R^^ is the Ricci tensor, and *R{AY^^ is the Hodge dual of the field strength of the U{1) 
transformations. Note that the T^-terms in TZ{M)^^^'^ cancel exactly the contribution from 
f^'^. Finally, Lcm denotes additional terms belonging to a compensating multiplet, that add some 
missing degrees of freedom needed for Poincare supergravity, and take care of the D field ap- 
pearing linearly in the Lagrangian. We discuss two equivalent possibilities: a nonlinear multiplet 
and a hypermultiplet. 

The bosonic content of the nonlinear multiplet is a real vector field VJj, a complex antisym- 
metric SU{2) triplet of scalars Mij, and an SU{2) matrix of scalars These are subject to 
the constraint: 

-^R-D + VVa - ^VVa - \Mi,M'^ + = . (6) 

In the absence of i?^-terms this constraint may be imposed by 

L,^ = i{X'Fi - X'Fj) (^-^R-D + P'^K - ^V^Va - \m,,M'^ + D-^^^^D.^^) . (7) 

This cancels the D field in the first line of the Lagrangian ([2]), and leaves —i{X^Fi — X^Fi)R/2 
as the Einstein-Hilbert term. In the presence of i?^-terms the above Lcm does not suffice, and 



^We have assumed the /C-gauge fixing which wiU be defined later (jTU]). 
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one must also substitute the D field appearing in the hatted fields ([3]) in the Lagrangian using 

Alternatively, one may use a hypermultiplet as the compensating multiplet [6]. The bosonic 
content of the hypermultiplet are scalars with F, A = 1 . . . 2r and r is an integer. We have 

Lcm = -hlr^VaA^V^A'^ + ]^x{\r + D) , (8) 

where Ci^a is a covariantly constant antisymmetric tensor, and the hyper-Kahler potential x is 
defined by 

X = I^vaA.'A^'' ■ (9) 

In order to obtain Poincare supergravity one gauge fixes the bosonic fields (in addition there 
is a gauge fixing of fermionic fields) : 

i^-gauge : ba = 
A-gauge : X° = X° > 

gauge : ^\ = 5^ (nonlinear multiplet) 

A"^ = ^r^pT^ (hypermultiplet) , (10) 

where ba is the connection of the dilatations. Note that we have left out the Z)-gauge, which will 
be discussed later. 



3 Large Charge N = 2 Black Hole Solutions 
with i?^- Terms 

We assume a prepotential of the form: 

where Dabc, Da are constants, Dabc is symmetric in all indices, and A,B,C = 1 . . . Ny- The 
second term describes couplings in the Lagrangian. It may arise as a Qs correction in the 
large volume limi1§| of type IIA string theory compactified on a Calabi-Yau three-fold, or as an a' 
correction in heterotic string theory compactified on K3 x T^. We will treat the higher curvature 
terms in the Lagrangian as a small perturbation, in the spirit of This is valid for the exterior 



^The large Calabi-Yau volume approximation requires Im(X ^ 1. 
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region of black hole solutions in the large charge approximation. The physical dimensionless 
expansion parameter is one over charge squared: Q~^. It is however convenient to express this 
as an expansion in \Da\ ^ \Dabc\ (see also [3 HI [5]). To make this explicit, we have inserted 
the expansion parameter e in front of the second term in ( fTTj) . and at the end we will set e = 1. 

We look for static spherically symmetric solutions, where the metric takes the form: 

ds^ = _e-2t^M f[r)de + e^^^'-) [firY^dr^ + sin^ Odcj)^ + r^dO^) . (12) 

Introduce the rescaled variables: 

Y' = e^X^ 
T = e^^I 

e"^ = i{Y^Fi{Y,T)~Y^Fi{Y,f)) . (13) 

The latter is called the Kahler potentiallfl 

Consider black holes with one electric charge go and p"^ {A = 1,2,3) magnetic charges. In 
our convention V = D^bcP^P^P^ > oE For go > one can have a supersymmetric solution 
[ill [21 [9]. By reversing the sign of the charge go, but taking the moduli to depend on the absolute 
value, one can have an extremal non-supersymmetric solution [121 IH E]- At the i?-level, i.e. 
without i?^-terms, in addition to the sign changes, the supersymmetric and non-supersymmetric 
solutions differ also in the form of the auxiliary T field. The thermodynamic properties of the 
non-supersymmetric solution can be obtained by an analytic continuation. This is no longer true 
when including higher curvature corrections. We will construct the i?^-level solutions for both 
cases. In the following we will denote E = D^p^. 

We are interested in black hole solutions of the curvature corrected theory, to first order 
in e. As a starting point for our ansatz, we may take the i?-level solution (e = 0), with the 
prepotential F(e = 0) replaced by F{e). This, however, proves to be insufficient, and we need to 

^The scaling X{zy = f-'^^'^X^ used in some previous works is not general enough. 
"^The common convention uses D < and a reversed sign for go- 
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introduce a further general linear e-correction to the fields. We look for solutions in the form: 

— —-y^{l-\-e^A{T)) (no summation) 



2 y uq 

where {ko, k^) > are constants with either k^ = k^ = k^ or otherwise only 1^333 7^ 0, and 

yo - ^ + ^- (15) 

We use 

T.-4e-ra--4(-flj+sg.ta)^)^i + 0(e) , (16) 

which is a sufficient approximation since T always comes with a factor of e. 

The event horizon is located at r = 0. In order for the perturbative expansion to be valid, 
we require |e^(r)| <^ 1 for all ,^-functions for r > 0. In addition we set the boundary conditions: 
limr^ooC(^) = 0- This gives an asymptotically flat solution. 

Let us introduce the dual field strength: 





+ — 
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GaM = 2i^^ = FijF;,' + ^{Fj- FuX')T;, + F^,F;, , (17) 
where we have considered only bosonic terms. Due to spherical symmetry we have 

^231 = -^^017 ■ (18) 
The field strengths F^/ may be extracted from the following equations: 

2(ImF,j)Fo-/ = 6-23/ - FjjF^, + llm( (F, + FjjX' - 64F^,(2Coioi - D)) Tol) , (19) 

where we used spherical symmetry and 

nMr%, = C-"^, + DS^^S^l - 2S\;*R{Ar\^ , (20) 
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where C"''^^ is the Weyl tensor. The magnetic parts of the field strengths are obtained from 
Bianchi identities, which for our static spherically symmetric metric f|T2|) give: 

^23 - P 

G231 = ^e-^V. (21) 

For our choice of charges, and the complex-valued form of the prepotential (fTT!) and the ansatz 
(fill we get 

Fo' = 7^ 0^230 - ^FoAFi + ^ (Fq + FojX^ - 64F;fo(2Coioi - D)) T,, 
For the U{1) and SU{2) connections we assume 

Aa = 

Vy = 0. (23) 

The equation of motion for the SU(2) connection is always satisfied by the vanishing SU{2) 
connection, for a bosonic background and with either the y-gauge for the nonlinear multiplet 
(1 101) or covariantly constant hypermultiplet scalars. This is because the SU{2) connection and 
its derivatives, appear then in the Lagrangian ([2]) always in at least a quadratic form|§ The 
vanishing of the SU{2) connection implies also Y^^ = ^13] . 

When using the nonlinear multiplet, the auxiliary field D may be determined by the constraint 
on the nonlinear multiplet ([6]). We assume 

Va = 

M,, = 

^\ = (24) 

where the later equation is the V^-gauge. The equations of motion for Mij and are trivially 
satisfied by the above assumption, where for the latter we assume a vanishing SU (2) connection. 
We therefore remain with the constraint: 

D = -IR- (25) 



^More generally, the SU{2) field equations arc automatically satisfied since the solution is a singlet under the 
SU (2) symmetry. 
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When using the hypermultiplet, we assume covariantly constant hypermultiplet scalars: 



VaAl = . (26) 

The equation of motion for the hyperscalars then gives (l25l) . Solving the other equations of 
motion we will find that for all our cases: 

X = -2 + 0(e2). (27) 

For our ansatz to constitute a solution, it must satisfy the equations of motion for the metrical, 
the moduli F^, the auxiliary field Tq]^, the U{1) connection Aq, and either Va for the nonlinear 
multiplet or the auxiliary field D for the hypermultiplet. In the cases that we solved, we observed 
that when using the hypermultiplet, the equation of motion for D has an overall factor of {k^ — 
koY, after substituting the ansatz. For = ko, one has to take this limit only after solving the 
equations of motion, in order not to lose a constraint. 

From the Einstein-Hilbert term in the Lagrangian ([2]), one sees that "Newton's constant" is 
given by the unsealed Kahler potential: 

G],' = e-^ = i (x'Fj{X, A) - X'FiiX, A)) = 1 - e^u{r) . (28) 

Usually one fixes = 1 a.s the dilatational D-gauge choice. This, however, is too restrictive and 
does not always allow a solution. Therefore Gjv is a function of the radial coordinate, resembling 
the case of dilaton gravity. The metric in the Einstein frame is given by 

9^ = G-^^g^, . (29) 
The ADM mass (in Planck units) for a non-normalized metric is given by 

^.f I™ = 9fA^) (l - 9fr{^r"^-f- + 0(^)) . (30) 

One may see this by applying the coordinate transformation t {—gft (oo))~^''^t, r — > gf^.{oo)~^^'^r^ 
to get the conventionally normalized line element. 

The central charge is the conserved charge associated with the graviphoton: 

Z = — (f e^'\F-^Fi{X,A) - G7X^) = — I e^'^{F^ Fi{X,A) - GiX^) = 

= lim e^/2(/F7(X, A) - qiX^) = lim e^/2(/F/(F, T) - qiY^) , (31) 



^For a discussion on the derivation of the metric field equations, see [T4i appendix B]. 
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where the first equahty is valid due to the N = 2 supersymmetry of fiat space at infinity. This 
also allows to express Z in terms of the auxiliary T fielcu: 

Z = -— / e-'^/'T- = -\ hm e-'^WT- . (32) 
IGvr Jg2 4 r^oo 

Z is determined by the charges and the asymptotic moduli values at infinity, and does not receive 
higher curvature corrections. It reads 



\Z\ = ^{\ho\DABch^h''h''y/^\3e + sgn(go)A;o| = ]\3k^ + sgn{q,)ko\ 



(33) 



where ^ 

{hQ,h^)= (no summation) , (34) 

and where in the second equality of (!33|) we imposed the normalization gtt{oo) = —1: 

MholDABC'h'^h^h^ = 1 • (35) 

The supersymmetry algebra requires M > \Z\. 

Our calculations were done using Maple with GRTensor. 



3.1 Supersymmetric Black Holes 

In the supersymmetric case (go > 0) with = k"^ = k^, the solution reads 

8Er{k^f{k^ - ko){9rk^ - rko + Sk^ko) 



V{r + k^)^{r + ko)^ 



4:Er{ey{e - A;o)P(r) 
^^^""^ ~ V{r + k^Y{r + koY{?,rk^ + rk^ + Ak^k^) ' 



where 



P(r) = 144r^A;^ - 16r^A;o - 27r'^{k^f + 22Qr^k^ko - Tr^fcg - 100r{kYko + WOrk^k^ - 6A{k^)'^kl . 

(37) 

This can be obtained either from the full second order equations of motion, or from the first 
order = 1 supersymmetry equations [6j similarly to [oj^. At the horizon (r = 0) we have: 

^u{0) = ^t{0) = . (38) 



^The T field differs from the graviplioton e'^{F ^ Fj — Gj X^), when i?^-ternis are included. 

^The case discussed in [5] with Dabc — -D123 and constant y^, F^, is not soluble in our approximation scheme. 
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128-E \ 

S = 27T^qoV ( 1 + — — + 0(Q-2) ^ 27r^qo{V + 256E) . (40) 



The horizon solution agrees with previous results [2]. 

The entropy is given by the Wald formula for the supersymmetric case with i?^-terms [2] : 

^ = J'i'o(4^"^^-^^(^^l^o"^l')) ' ^^^^ 
where A is the area of the event horizon. Plugging in the solution gives the expected result: 

In the supersymmetric case the ADM mass saturates the BPS bound: 

M = \Z\ = ^{Se + ko) . (41) 

This result is exact and is identical to that of the i?-level. 

Note that in the special case = ko, the ^-functions vanish in all space. This means that 
the i?^-level solution is simply the i?-level ansatz with F{e = 0) replaced by F{e). One may ask 
whether this behavior continues to higher orders in e. Assuming that the prepotential (|TT|) itself 
does not contain higher orders in e, the solution to order would be 

Q5536E^rkl 
^ 32768EVfcg(7r-fco) 



3.2 Non-supersymmetric Black Holes 

For simplicity, we consider the cases with Dabc = -D333 or Dabc = -D123, and Da = -D3 and 
k^ = k^ = k^ = ko. At the i?-level, the last condition gives the non-supersymmetric version of 
the double-extremal black hole. We were able to find similar solutions for different combinations 
of Dabc^ and Z^a's, providing that for each term such as there is a at least one corresponding 
term Dab3 (this was not required in the supersymmetric case). We then give the generalization 
of the Dabc = -D333 case to arbitrary k's. 
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In the non-supersymmetric case (go < 0) with — — — ko, the solution reads 



V{r + ko)^ 
SEkoi'ir^ + 25r^ko + 60rk^ + 30k^) 
15V{r + koY 

AE (-120r(r + A;o)=^ In ^ + 118r=^A;o + 2^1r'^kl + 22Arkl + 240A;^) 

lbV{r + k^Y 

AE (I20r(r + A;o)=^ In ^ - lld>r^kQ - 297r^k^ + 736rA;g - 120k^) 



15V{r + koY 

For Dabc = -D333 and Da = D^: 

4£; (-120r(r + kof In ^ + lOSr^ko + 267r^kl + 19Ark^ + 27Qk^) 
" A5V{r + koY ' 

and ^i(r), ^2(^) are irrelevant. For Dabc — -D123 and Da — D3: 

AE (-840r(r + kof In ^ + 856r3A;o + 2169r'^k^ + 1658rkl + ISO/c^) 



15V{r + koY 

AE (360r(r + kof In ^ - 37Ar^ko - QSlr^fc^ - 732rA;g + 60A;^) 



At the horizon (r = 0) we have, 



For Dabc = -D333 and Da^ D3: 
For Dabc = -D123 and Da^ D3: 



m -- 
m -- 

&(o) -- 

m 



15V{r + kof 



6AE 



V 

16E 

~V~ 
ME 

32E 
V~ 

2AE 



V 



6(0) = ^ 

6(0) = 6(0) = ^ 
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The latter case agrees with [5], where the horizon solutions were obtained by extremizing the 
entropy function^ Note that the radial derivatives of our solutions for ^Air),^oi^)^^TiT") diverge 
at r = 0. The curvature, however, is regular. 

The entropy formula that we use in the supersymmetric case fl39l) . is no longer valid in the 
non-supersymmetric case. In the latter case, one can use either Sen's entropy function method 
[U [5] , or the Wald entropy formula for the non-extremal case with i?^-terms [H] : 



S = lim - AA . Im(F^(|Tol|^ + I6C0101 + 16D))) 



(49) 



In fact, in our solutions Coioi,-D do not contribute to the entropy to first order in e, and thus 
the formula does reduce to the supersymmetric one. Note also that GAr(O) 7^ 1. Plugging in the 
solution gives: 

40E 



S = 2n^-qoV (^1 + + 0{Q~^) , (50) 

as in [U E]. As discussed therein, this differs from the statistical microscopic entropy [15], due 
to higher curvature D-terms, which are not taken into account here. Using the entropy function 
method with our first order solutions, one can get the entropy to second order: 



40E 8576E 



S = 27ry/-qoV + — ^ J + 0{Q-') . (51) 

The central charge reads 

\Z\ = ^{-qoV f = ho , (52) 



and the ADM mass takes the form: 



(53) 



Here the mass gets an curvature correction. This correction comes from the functions 
Ol'")' ^o('"); ^a('"), which vanished in the supersymmetric case. The higher curvature correction 
changes the mass in an opposite direction to the entropy. For the mass to decrease as conjectured 
in [16], V and E must have the same sign. This would also mean that the entropy increases. 

In the non-supersymmetric case (go < 0) with Dabc = -D333 and Da = -D3 and arbitrary fc's, 
the solution is given in the appendix. The horizon limit, and the entropy are the same as in the 
= ko case. The central charge reads 

|Z|— I 

I I ^1 



\3k^ - ko\ , (54) 



^Note that [5, uses a different scaling than 
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and the ADM mass takes the form: 

1 AEk^ko ( 12(efko In f - (e - ko) (3(k^f + Ui^fko - hk^'kl + kl)] 

M = l(3.'+.o)+ ^ ^^f^, Uo(g-^) . 

(88) 

The mass correction term has the same behavior as in the k^ — ko case. 
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A The Non-supersymmetric Solution with Arbitrary /c's 

In the non-supersymmetric case {qq < 0) with Dabc = -D333 and Da = -D3 and arbitrary /c's, the 
solution reads 



8Ds{k^f{9r\k^f - 2r^k^ko + r^fcg + 16r{k^)^ko + 8(^3)2^2^ 

A333(p=^)2(r + fc3)4(r + A;o)2 
leDsik'fko (2k'{r + k^fP,{r) In g|±^ + ^(^3 _ ko)P2{r)) 

i^333(p')V(r + A;3)2(A;3-A;o)5 
8£>3(-2(r + k^)\r + h)L^{r) + r{k^fko{k^ - ko)P,{r)) 

3A«3(p^)V2A;o(r + k^f{r + ko){k^ - kof 
4D3(2(r + e)^Lo{r) + rkH-^k^ ko)Ps{r)) 
D333{p^)^r^k^ko{r + k^Yir + ko){k^ - kof 



^ . , 4A3(2(r + efLrir) + rk^koik^ - ko)Pi2{r)) 

^^^^^ Ds33{p^yr^k^ko{r + k^y{r + koy{k^-ko)'^{3rk^-rko + 2k^ko) ' ^ ' 

and ^\{r)^^iif) are irrelevant, and where 

L3(r) = {k^klP^ir) In f f + + T\k^fP,{T) In + T\k^ + k^W - In 

- r3(A;3 + ko){k^ - ko)%rk^ + rifco + 2k^ko) In 
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Lt{t) = 2{e)X{r + ko)P,{T)\n y 

ko{r + k-^) r + ko 

T —I— 

+ 2r\r + ko){k^ + ko){k^ - hfPii{r) In 
Pi(r) = 2r{k^f + Srk^kQ + 2rkl + ?,{k^fko + ^k^kl 

P^{r) = ~12r\ey - 12r^eko - 20r(P)3 - 25r{efko - Ark^kl + rkl - Q{k^f - 18{k^fko 
P^[r) = -6r^{kY - 6r^{k^fko + 66r^k^kl + 6r^kl + 3r{k^fko + SOrik^k^ + 27rk^kl + 

+ {e)% + mkyki + 9{efki 

P^(r) = -{kY + 5{k^fko + 9{k^fkl-21k^kl-12k^ 

P^{r) = 2r\k''Y - 8r\k^)% + 32r\k^)X + QOr^k^k^ - 6rX + 4:r\k^)^ - 2lr\k^)% + 
+ 103r^(A;=^)^A;^ + U^r^k^kl + hlr^'k^kl - 6r^A;^ + 2r\k''f - l2r\k^fkQ + 
+ %%r''{k^fkl + 162r2(A'r^)'^A';jj + l2&T''{kykl - 6rVk^ + Wr{k^)%o - 26r{kYkl + 
+ 189r{kYkl - 9r{efk^ + Abri^fk^^ - 9rekl + 2{k^fkl + 2Q{k^fkl + IS^k^kl 

Peir) = 2r^{k^f - 12r^{k^fko - ASr^k^k^ - 2r^kl - r^k^ - 15r\k^fko - 42r\k^fkl + 

- 59r2pA;3 - 3r^k^ - 2r{k^fkQ - 22r{k^fkl - 38r{k^)^kl - ISrk^k^ - {k^fkl + 

- 10(A;3)3A;3 - 9{kykt, 

P^{r) = r{ky - Ar{k^fko - 2r{kYkl + 36rk^k^ + 9rk^ + 2{kYko + 2{k^fkl + 6{k^fk^ + 
+ ISk^ko + 12kl 

Psir) = 2r\k^)'' - 12r\k^)% + mr^k^kl + IQr^k^'fkl - I9>r\k''fkt, + d>r^k^kl - 2r% + 
+ 4r3(P)7 - 2br^{k'^fko + ZQbr^^k'^fkl + 22Sr'^{k'^fkl + AQr^^k^k^ + Ur'^^k^fkl + 

- hr^k^k^ + 2r2(P)8 - Qr'^{kykQ + 2Ur'^{k'^fkl + 'HQr'^ {k'^f kl + USr'^ {k'^)% + 
+ m-^{k'^fkl - 4r2(P)2A;6 + 22r(A;3)8A;o - Qrik'^fkl + 354r(A;3)6A;3 - 34r(A;3)5A;4 + 
+ 80r(P)^A;^ - 16r(P)3A;^ + 4(Ar^)*^A';^ + 40(A;2)^A;ij + 36(A;^)*^A;;5 

P^{r) = -6/(A;=^)^ - 4r^(A;3)3A;o + 9Qr\efkl + 36r^A;3A;3 - 2r^A;^ - br^efko + Slr\efkl + 
+ 249r^(A;^)2A;o^ + 35r3A;^A;^ + Zr^k^'fk^ + 46r2(A;3)4A;2 + lUr^k^'fkl + 182r2(A;^)2A;4 + 
+ 15r^A;^A;^ + hr{k^fkl + 55r(A;^)^A;^ + 95r(A;=^)^A;^ + 4hr{kykl + 2{k^fkl + 
+ 20(A;=^)^A;^ + I9>{k^fkl 

P^^(r) = -r2(A;3)5 + 5r\e)% + 8r2(A:3)3A;2 - 28r2(A;3)2A;3 - 27r^k^k^ + Sr^k^ - 2r(A:3)5A;o + 
+ 8r(A;3)^A;o' + iri^fkl - 72r{k^fk^o - ISrk^kl - 2{k''fkl - 2{kYkl - Gi^fk^o + 

- ISik^k^o - 12k^k^o 

P^^{r) = r\kY + r^kl + 2r{kYko + 2rk^kl + 2{kYkl 
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Pi2(r) = -Ar''{ky + 16r^(A;=^)^A;o - b^r^k^fkl - 2&^r\k^Ykl - 2Qr''{k^fkt, + 2Ar''{k^fkl + 

- 16r^A;=^A;^ + Ar'kl - I2r\k^f + 50r^(A;=^)^A;o - 214r^(A;=^)^A;^ - I22&r\k^fkl + 

- 724r^(A;3)^fc^ + mr^{kykl - ^r^ik^fkl + 2r^k^kl + Ar^k^ + lZ2r^{k^f + 

- 678r^{k^fko + 1166r^(fc^)^fc^ - 3778r\kYkl - 2062r^{k^)^k^ - A50r^{k^)^k^ + 

- 86r\kyk^o - 22r\kykJ, + ISr^k^k^Q - 31r\k^y° + 586r^(A;^)^A;o - 2706r\k^fkl + 
+ 26Ur\k^ykl - 9328r\kyk^ + 822r\k^fkl - 726r^(P)^A;^ - 70r^(P)^A;^ + 

+ 39r\k^fkl - 70r^{kY°ko + 678r%kykl - 3926r%k^fkl + UbOr^k^k^ + 

- 73Mr%k'fk'o + 1922r'{k'fk^o - QlOr^ik^kl + BOr^k^k^a - SOr^k^^l + 
+ 176r\k^)X - 2SA2r\k^)X + 152r\kYkl - 232Ar\k^fkl + 728r\k^fkl + 

- l3Qr\k^Ykl - AQrik^'fX - 2mr{kykl - 776r{e)X - 200r{kYk^o + 

- 80r{kykl + 16r{kYk^o " 8(A;=')'°/t^ - 80(^)^1 - 72{k^)X ■ (57) 
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